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1. INTRODUCTION 

Edge detection is one of the fundamental tools to find the boundaries of an objects for object 
segmentation and feature extraction [1]. There are two main edge detection operators used for edge detection; 
the first derivative based and the second derivative based. The first derivative edge detector operator 
computes the image gradient values to detect image edges. The edges are estimated based on gradient 
magnitude which is calculated in the x and y directions [2]. The second derivative edge detector operator uses 
zero-crossing for edge detector. The method finds the rate of change in grey intensity and then detects the 
local maxima in the gradient magnitude. It locates the centres of thick edges and the localization of edges is 
good [2], [3]. 

In most cases, the edge pixels are discontinued in the edge detection. This is because the process is 
based on local change in image intensity. In practical, a connected curve that shows the boundary of an object 
is desirable for feature extraction. Generally, this type of edge detector is sensitive to noise, and this is 
because the differential edge detectors behaves as a high-pass filter and thus, it has a tendency to amplify 
noise [4], [5]. Due to this, it is difficult to differentiate the true and false edges based on the differential edge 
detector. Multiscale edge detection has attracted a lot of attention in image processing applications due to the 
fact that edges are multiscale in nature [6]-[8]. The edges extracted from a different scale or resolution are 
then combined to get the final edge detection result. Multiscale edge detection method has the ability to 
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describe the variety of the edge structures [8]-[10]. Wavelet transform is one of the methods often used for 
multiscale edge extraction. The wavelet coefficients represent the scale of features in an image. It can be used 
to measure how closely the correlated wavelet is with each section of the image [11 ]-[13]. 

S. V Seeri et al. used Haar discrete wavelet transform as multi-resolution method to model the 
characteristics of textured images [14]. This method well characterized the edges of the textured images. The 
method fused edges that were extracted from three sub bands which were horizontal (H), vertical (V) and 
diagonal (D) for edge information. After the edges were extracted, the method also integrated with Sobel 
edge detector, fuzzy thresholding and morphological operators method to segment and classify the text 
regions. One of the advantages of the method is that it is able to localize the text regions accurately within the 
three types of scripts which are English, Kannada and Hindi. Wavelet transform, multiscale method is 
sensitive to local regularity and it has the limitation of detecting orientation of singularity curve. K. Guo et al. 
applied framework of Shearlet transform into wavelet transform to improve the edge orientation and curve 
point detection [15]. The modified method had the ability to estimate the local flatness of an edge and it 
could also provide detailed information on geometry edge. Obviously, the method has the advantage of non- 
parabolic scaling and it is able to discriminate features. 

Thigh bone is the longest and strongest bone of the entire human body. Various diseases often occur 
in children and adults. Magnetic resonance imaging (MRI) is an important imaging tool used to investigate 
the presence of certain abnormalities in thigh [16], [17]. Because of the nature of quantitative MRI imaging 
data, this requires a considerable amount of human intervention and expertise training opportunities. This 
process is tedious, time consuming and labour intensive. The development of automatic computer-aided 
diagnosis system has the additional benefits to help in diagnosing various pathologies [18], [19]. The purpose 
of this study is to present the formation of multiscale edge detection method that uses three different kernel 
sizes of averaging filter to smooth image into three different scales. The edges were extracted from three 
averaging filters that represented coarse to fine scale. Three classical edge extraction operators namely 
Prewitt, Sobel and Laplacian were used to extract the edges to establish the multiscale edge detection 
method. The performance of the three classical edge extraction methods was compared. 

This paper is organized as follows: in section 2, theories related to averaging filter and edge 
detection methods on first-order and second-order edge detector operator are presented, while in section 3, an 
edge extraction method for multiscale edge extraction method is proposed. To clarify, we used three edge 
extraction operators for the process of multiscale. The experimental results and discussions on the 
comparison results are then discussed in section 4. Finally, some conclusions are drawn and discussed in 
section 5. 


2. AVERAGING FILTER AND EDGE DETECTION METHOD 

Generally, averaging filter is also known as mean filter. It has smoothing effect by average value of 
the image to reduce noise. On the other hand, classical edge detector is often used for edge detection due to 
its simplicity [20], [21]. 


2.1. Averaging filter 

The implementation of average filtering is to smooth images. The average filter reduces the amount 
of intensity variation by averaging the pixel itself with the pixel value around its neighbours [22], [23]. 
Figure 1 shows the kernel averaging operator for 3x3, 5x5 and 7x7. A 3x3 square kernel as shown in 
Figure | (a) is often used to smooth images. The smoothing effect produced depends on the size of the 
averaging filter. Some details in the image will be smoothed by using a 3x3 averaging filter. The effect of 
smoothness will be increased by using a 5x5 averaging filter as shown in Figure 1 (b). Next, Figure 1 (c) 
show the kernel of 7x7 averaging filter. The smoothing output of a 7x7 averaging filter will be more than a 
5x5 averaging filter. Therefore, larger kernels can be used for more severe smoothing. 


(a) (b) (c) 





Figure 1. Averaging filter, (a) 3x3 kernel, (b) 5x5 kernel, (c) 7x7 kernel 
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2.2. The first order edge detector operator 

Figure 2 shows the Prewitt edge operator model [24], Figure 2 (a) shows vertical direcation, 
Figure 2 (b) shows horizontal direction. Prewitt edge has one by one model that is used to detect edges in an 
image. Based on the model, the values that are symmetrical and central (x, y) will make the model take the 
maximum value that is similar to the detected region. 

Figure 3 shows the Sobel edge operator model [25], Figure 3 (a) shows vertical direcation, 
Figure 3 (b) shows horizontal direction. Sobel operator computes an approximation of edge from the gradient 
of the image intensity function. Sobel operator gives weight to the point that lies closer to (x, y). As a result, 
it has smoothing effect and it is also less sensitive to image noise. 





Figure 2. Prewitt edge model, (a) vertical Figure 3. Sobel edge model, (a) vertical direcation, 
direcation, (b) horizontal direction (b) horizontal direction 


2.3. The second order edge detector operator 

Laplacian is the second spatial derivative that has a 2D and isotropic measure of an image. Figure 4 
shows the two commonly used Laplacian edge operator models [25], Figure 4 (a) shows outward edges, 
Figure 4 (b) shows inward edges. The operator enhances the regions of intensity discontinuities and 
highlights the regions that have rapid intensity change. Thus it is often used to detect fine edges. Due to its 
characteristic that it is sensitive to noise, the image often is smoothed with Gaussian smoothing filter to 
reduce noise. 





(a) (b) 


Figure 4. Laplacian operator, (a) outward edges, (b) inward edges 


3. THE PROPOSED MULTISCALE METHOD 

The method adapted the concept of multiscale method to extract edges from different resolutions on 
thigh images. The proposed method is illustrated in Figure 5. Three different kernel sizes of averaging filters 
were used to smooth the sharpened images and produce three different scales. This was conducted in order to 
achieve the coarse to fine scale for the multiscale concept. The method involved three main processes which 
were image sharpening, multiscale and edge extraction by classical edge detectors as shown in Figure 5. 

Firstly, the thigh images were sharpened by using Prewitt operator which rotated in eight directions; 
north, west, south, east, north west, south west, south east and north east. In this process, eight sharpened 
images were produced. Prewitt operator was chosen due to its its less computation complexity. Prewitt edge 
detector produces edges which were the most similar to the original images because of its weight on the 
model. 

The concept of multiscale method which included coarse-to-fine edge extraction method based on 
the size of the averaging filter was then applied to smooth the sharpened images. The eight sharpened images 
were then smoothed with three different averaging filters which were 3x3, 5x5 and 7x7. The averaging filter 
was used to produce images from three different scales which were high scale (3x3 kernel), middle scale (5x5 
kernel) and low scale (7x7 kernel). Edge extraction was then performed on the output of each averaging filter 
using three different classical edge extraction detectors. The edges extracted from the eight output images 
from each averaging filter (Prewitt, Sobel and Laplacian) were then combined. A comparison analysis was 
carried out based on the edge detection from each edge detector. Figure 6 illustrates the steps of the edge 
extraction and the edges combination for each scale. 

The comparison analysis was then carried out for each detector. Figure 6 shows that the output from 
the three averaging filters were channelled into Prewitt, Sobel and Laplacian edge detectors for edge 
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extraction from coarse-to-fine edge extraction. The output of multiscale edge extraction from each edge 
detector was then compared. 
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Figure 5. The concepts of multiscale edge extraction method 
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Figure 6. Comparison of multiscale methods 


4. RESULTS AND DISCUSSION 

Edge extraction is a process to detect the boundaries, so that objects can be extracted from the 
background image. Figure 7 (a) shows the original MRI image. The edges extracted from an edge detector 
operator often depends on image brightness whereby curved line or edge is extracted when the intensity 
variation is high in the image. In Figure 7 (b) and 7 (c), it is obvious that the edges extracted from Prewitt and 
Sobel edge detectors are very similar. The extracted edges from the outlines of the boundaries on muscle, fat 
and bone morrow are clear. However, it can be seen that some discontinued edges appear in the image. The 
discontinuity of surface depth and surface orientation cause discontinuity of edges. For Laplacian edge 
detection, the detector could extract more fine detailed edges as shown in Figure 7 (d). Unfortunately, more 
discontinued edges were found in the image and the fine detailed edges appear as noise in the image. 
Generally, the discontinuity of edges could fail to locate important edges that could cause inaccuracy of 
object extraction. 

The multiscale edge extraction method extracted edges over a wide range of scales. When the 
multiscale edge information fuse together, an edge map could be obtained. The edges information from 
coarse to fine scale that made the edges were robust. The proposed multiscale edge extraction method is to 
begin with image sharpening. Images were sharpened using Prewitt compass operator. This was because 
Prewitt was simple and the edges extracted were very similar to the original image. As shown in Figure 8 (a) 
to Figure 8 (f), eight images were produced based on eight different orientations. The sharpened images 
enhanced the edges of the thigh image from eight different orientations. 
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(b) (d) 


Figure 7. Edge detection for classical edge detector, (a) original, (b) prewitt extraction, (c) sobel extraction, 
(d) laplacian extraction 





Figure 8. Images sharpen on eight orientation using Prewitt compass operator, (a) east, (b) north, (c) north 
east, (d) north west, (e) south, (f) south east, (g) south west, (h) west 


Applying the multiscale concept, three differences sizes of average filters were used to smooth the 
sharpened images for three different scales. Figure 9 (a) shows the original MRI image. Figure 9 (b), 
Figure 9 (c) and Figure 9 (d) show the output of the averaging filter of 3x3, 5x5 and 7x7 respectively. The 
output image was a weighted sum of the input pixel. As the mask size increased, the image became blurry 
and it react to coarse-scale (low scale) methods. The smaller the mask size, it responded more to the spatial 
accuracy of fine-scale (high scale) methods. 





(b) (c) (d) 


Figure 9. Three difference of average filters, (a) original, (b) average 3x3, (c) average 5x5, (d) average 7x7 


Three different classical edge detectors were used to extract the edges from the eight orientations of 
the sharpened thigh images. The edges extracted from each sharpened orientation image from the same scale 
were then combined. Figure 10 shows the edge extraction comparison between Prewitt, Sobel and Laplacian 
edge extraction methods. Figure 10 (a), Figure 10 (b) and Figure 10 (c) show the output of Prewitt 3x3, 5x5 
and 7x7 respectively. Figure 10 (d) is the output of multiscale edge extraction for Prewitt edge detector. 
Figure 10 (e), Figure 10 (f) and Figure 10 (g) show the output of Sobel 3x3, 5x5 and 7x7 respectively. 


Comparative studies of multiscale edge detection using different edge ... (Belinda Chong Chiew Meng) 


19840 ISSN: 2302-9285 


Figure 10 (h) is the output of multiscale edge extraction for Sobel edge detector. Figure 10 (1), Figure 10 (j) 
and Figure 10 (k) show the output of Laplacian 3x3, 5x5 and 7x7 respectively. Figure 10 (l1) is the output of 
multiscale edge extraction for Laplacian edge detector. 





Figure 10. Result of three multiscale methods, (a) prewitt 3x3, (b) prewitt 5x5, (c) prewitt 7x7, (d) prewitt 
multiscale, (e) sobel 3x3, (f) sobel 5x5, (g) sobel 7x7, (h) sobel multiscale, (1) laplacian 3x3, (j) laplacian 
5x5, (k) laplacian 7x7, (1) laplacian multiscale 


The output of the multiscale edge extraction in Figure 10 (d) for Prewitt detector, Figure 10 (h) for 
Sobel detector and Figure 10 (1) for Laplacian detector could product continuous edges if compared to 
Figure 7 (b), Figure 7 (c) and Figure 7 (d). In Figure 10, the resultant multiscale edges can be seen 
continuous and the boundaries are clearly seen. In the experiment, both Prewitt and Sobel edge detectors 
produced the edges that were most similar to the original image. Prewitt operator was sensitive to horizontal 
and vertical edge; Sobel operator was sensitive to diagonal edge. Due to the proposal of multiscale method 
considering all eight directions for edge extraction, the combination of eight direction edges caused the 
Prewitt and Sobel operators to produce edges that were similar to the original image. 

Laplacian operator is the second derivative operator that the method searches zero crossing to detect 
edges in an image. As shown in Figure 10, when the kernel size gets bigger, Laplacian 7x7 (Figure 10 (k)) 
may yield better results in edge localization compared to Prewitt 7x7 (Figure 10 (c)) and Sobel 7x7 
(Figure 10 (g)). Besides this, in the comparison analysis, it could be seen that the Laplacian operator was 
very sensitive to noise as shown in Figure 10 (1). It was due to Laplacian operator attempting to find zero 
crossing in the process of extracting edges. The fine detailed edges were extracted and amplified that this 
caused the multiscale edge extracted become noisy by using Laplacian operator. The operator tended to 
accentuate noise and cause noises to be detected in the background. 

Figure 11 shows the comparison of other different parts of thigh. In the figure, it can be seen that the 
edge extraction from Prewitt and Sobel detectors was close to the original thigh image. However, Laplacian 
detection contained noise. This can be concluded that Prewitt and Sobel are more suitable for edge extraction 
on thigh image. The edge extraction preserves the similarity if compared to the original image. 
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Patient 1 


Patient 2 


Patient 3 





MRI image Prewitt Sobel Laplacian 


Figure 11. Comparison of multiscale for Prewitt, Sobel and Laplacian on other thigh images 


5. CONCLUSION 

This paper presents a multiscale method that was created by using three different kernel size of 
averaging filters. The averaging filter smoothed the image and produced three different scales according to 
the kernel size of the filter. Prewitt, Sobel and Laplacian classical edge detectors were then used to extract the 
edges at different scales. The performances of these three multiscale edge extractions were then compared. In 
this study, it can be concluded that Prewitt and Sobel operators could extract edges most similar to the 
original image. Laplacian could extract fine detailed edges. However, it is sensitive to noise. Laplacian 
operator could convolve with some smoothing filters to reduce noise in the background. Finally, it can be 
concluded that the proposed multiscale method that uses either Prewitt or Sobel edge detector as an edge 
detector offers a great advantage in term of producing multiscale edge detection. In addition, due to its 
simplicity, the method could reduce computation complexity and more importantly it has good edge 
preservation and could help in feature extraction. For further improvements, it is recommended that the 
image should be de-noised before the multiscale process takes place. 
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